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We report measurements of the lifetimes of the 3d ^05/2 and 3d ^^03/2 metastable states of a single 
laser-cooled *°Ca^ ion in a linear Paul trap. We introduce a new measurement technique based 
on high-efficiency quantum state detection after coherent excitation to the D5/2 state or incoherent 
shelving in the D3/2 state, and subsequent free, unperturbed spontaneous decay. The result for the 
natural lifetime of the D5/2 state of 1168(9) ms agrees excellently with the most precise published 
value. The lifetime of the D3/2 state is measured with a single ion for the first time and yields 
1176(11) ms which improves the statistical uncertainty of previous results by a factor of four. We 
compare these experimental lifetimes to high-precision ab initio all order calculations and find a 
very good agreement. These calculations represent an excellent test of high-precision atomic theory 
and will serve as a benchmark for the study of parity nonconservation in Ba^ which has similar 
atomic structure. 
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I. INTRODUCTION 



The lifetime of the metastable D-levels in Ca+ is of 
high relevance in various experimental fields such as 
optical frequency standards, quantum information and 
astronomy. Trapped ion optical frequency standards 
Q and optical clocks |2| are based on narrow absorp- 
tion lines in single laser-cooled ions. With transition 
linewidths in the 1 Hz range Q, Q- values (frequency of 
the absorption divided by its spectral width) of w 10^^ 
can be achieved. As the lifetimes of the D-levels in Ca+ 
are on the order of 1 s, yielding sub-Hz natural linewidths 
of the D-S quadrupole transitions, Ca"*" has long been 
proposed as a promising candidate for a trapped ion fre- 
quency standard Q. Such long lifetimes together with 
the ability to completely control the motional and elec- 
tronic degrees of freedom of a trapped ion Q make it 
ideally suited for storing and processing quantum infor- 
mation Q. In Ca^ a quantum bit (qubit) of information 
can be encoded within the coherent superposition of the 
Si/2 ground state and the metastable D5/2 excited state 
[3 with very long coherence times iSl . In astronomy, 
absorption lines of Ca+ ions are used to explore the kine- 
matics and structure of interstellar gas clouds [Tol [Tlj 
and the D-level lifetimes are required for interpretation 
of the spectroscopic data. On the other hand, in the- 
oretical atomic physics Ca^ is an excellent benchmark 
problem for atomic structure calculations owing to large 
higher-order correlation corrections and its similarity to 
Ba+. The size and distribution of the correlation correc- 
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tions make it ideal for the study of the accuracy of various 
implementations of the all-order method. The properties 
of Ba"*" are of interest due to studies of parity nonconser- 
vation in heavy atoms and corresponding atomic-physics 
tests of the Standard model of the electroweak interac- 
tion [13. 
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FIG. 1: Ca -level scheme with relevant transitions. 

Experimental investigations of long atomic lifetimes 
have profited enormously from the development of ion 
trap technology and laser spectroscopy. Early exper- 
iments on the measurement of the D-level lifetimes in 

4oca+ El Q El El 113 used large clouds of ions and 

the lifetime was determined by recovery of fluorescence 
on the UV-transitions (S1/2 - P1/2 or S1/2 - P3/2J see 
Fig. ^ after electron shelving in the D-states or by ob- 
serving UV fluorescence after driving transitions from 
the D-states to the P-states. Shelving in this context 
means that the electron for a certain time remains in a 
metastable atomic level which is not part of a driven fluo- 
rescence cycle. These lifetime measurements were limited 
by deshelving induced by collisions with other ions or the 
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buffer gas used for cooling. Similar results using the same 
techniques have been obtained in an ion storage ring |l8j| . 



Much more accurate results can be obtained by per- 
forming lifetime measurements with single trapped ions 
|i9, 20, 21, 22, 23] or strings of few trapped ions ^ and 
employing the quantum jump technique. This technique 
is based on monitoring the fluorescence on the S1/2 - P1/2 
dipole transition while at random times the ion is shelved 
to the metastable D5/2-state where the fluorescence falls 
to the background level. For observing fluorescence both 
the Si/2 - Pi/2 (397 nm) and the D3/2 - P1/2 (866 nm) 
transition have to be driven to prevent the ion from resid- 
ing in the metastable D3/2-state. Shelving to the D5/2- 
state is initiated by applying laser light at 850 nm (D3/2 
- P3/2) T?! or at 729 nm (S1/2 - D5/2) Statistical 
analysis of the fluorescence dark times yields the lifetime 
T. The most precise measurement using this technique 
was carried out by Barton et al. ^2211 who found the result 
of r=1168(7) ms. 
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Here, we introduce a measurement technique |23 based 
on deterministic coherent excitation to the D5/2 state 
or incoherent shelving in the D3/2 state, followed by 
a waiting period with free spontaneous decay and fi- 
nally a measurement of the remaining excitation by high- 
efficiency quantum state detection. During the waiting 
time all lasers are shut off and no light interacts with 
the ion. This method basically is an improved version 
of a technique that was used earlier to measure the D3/2 
metastable level lifetime in single Ba"*" ions [2BI. The 
main advantage of this " state detection" method is that 
no residual light is present during the measurement which 
could affect the free decay of the atom. In addition, it al- 
lows for the measurement of the D3/2 level lifetime which 
otherwise is inaccessible with the quantum jump tech- 
nique. There exist only a few reported D3/2-level lifetime 
results for Calcium [1511161 fig but none from a single ion 
experiment. 



Figs. I^and H] compare the different exp erimental 

IIllllImHOplSp^Ell^ and the- 

oretical [23, 128L I29I I30I I311 13^ methods and results for 
the D5/2- and D3/2-level lifetimes. From Fig. [21 it is 
evident that the single ion measurements are the most 
accurate ones. Generally, lifetime measurements on sin- 
gle ions or crystallized ion strings are more accurate as 
systematic errors, e.g. due to collisions, can be reduced to 
the highest possible extend. Therefore, single ion D-level 
lifetime measurements for Calcium are of special inter- 
est. The existence of accurate D-state lifetime values is 
of special interest for theory as well since most studies 
of alkali-metal atoms were focused on the measurements 
of the lowest nP-state lifetimes and D-states are much 
less studied. The properties of D-states are also gener- 
ally more complicated to accurately calculate owing to 
large correlation corrections. 



FIG. 2: Theoretical and experimental results for the D5/2- 
level lifetime. 
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FIG. 3: Theoretical and experimental results for the D3/2- 
level lifetime. 



II. EXPERIMENTAL SETUP AND METHODS 

For the experiments, a single ''^"Ca"'" ion is stored in a 
linear Paul trap in an ultra high vacuum (UHV) envi- 
ronment (10~^^ mbar range). The Paul trap is designed 
with four blades separated by 2 mm for radial confine- 
ment and two tips separated by 5 mm for axial confine- 
ment. Under typical operating conditions we observe 
radial and axial motional frequencies (o^ri, 0^1.2, i^ax) = 
27r(4.9,4.9, 1.7) MHz. 4°Ca+ ions are loaded into the 
trap using a 2-step photoionization procedure [13 • The 
trapped ^"Ca"*" ion has a single valence electron and no 
hyperfine structure (see Fig. . Doppler-cooling on the 
Si/2— Pi/2 transition at 397 nm puts the ion in the Lamb- 
Dicke regime [^ [^ . Diode lasers at 866 nm and 854 nm 
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prevent optical pumping into the D states during cool- 
ing and state preparation. For coherent excitation to the 
D5/2 state we drive the S1/2 to D5/2 quadrupole transi- 
tion at 729 nm. A constant magnetic field of 3 G splits 
the 10 Zeeman components of the S1/2 ^ ^5/2 multiplet. 
We detect whether a transition to D5/2 has occurred by 
applying the laser beams at 397 nm and 866 nm and 
monitoring the fluorescence of the ion on a photomulti- 
plier (PMT), i.e. using the electron shelving technique 
|34| . The internal state of the ion is discriminated with 
an efficiency close to 100% within approximately 3 ms 
(3^. The following stabilized laser sources are used in 
the experiment: two frequency-stabilized diode lasers at 
866 nm and 854 nm with linewidths of ~ 10 kHz and 
two Ti:Sa lasers at 729 nm (< 100 Hz linewidth) and 
794 nm (< 100 kHz linewidth), of which the 794 nm laser 
is externally frequency doubled to obtain 397 nm. The 
experimental setup and the laser sources are described in 
more detail elsewhere (2i, 1^ ■ 

III. MEASUREMENT OF THE D5/2 STATE 



LIFETIME 



A. Measurement procedure and results 
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FIG. 4: Simplified pulse scheme for the D5/2 lifetime mea- 
surement: the preparation consists of Doppler cooling, re- 
pumping, and optical pumping (2ms); followed by a 7r-pulse 
(few /is) and a detection periods (3.5 ms). The waiting time 
is varied between 25 ms and 5 s. 



The measurements consist of a repetition of laser pulse 
sequences applied to the ion. The sequence generally is 
composed of three steps (see Fig. 0J): 

1. State preparation and Doppler cooling, consisting 
of 2 ms of Doppler cooling (397 nm and 866 nm light), 
rcpumping from the D5/2 level (854 nm light) and op- 
tical pumping into the Si/2(ni=-l/2) Zeeman sublevel 
(397 nm polarized light). 

2. Coherent excitation at 729 nm with pulse length and 
intensity chosen to obtain near unity excitation (7r-pulse) 
to the D5/2(to = —5/2) Zeeman level. 

3. State detection for 3.5 ms by recording the fluores- 
cence on the Si/2 - P1/2 transition with a photomulti- 
plier. Discrimination between S and D state is achieved 
by comparing the fluorescence count rate with a thresh- 
old value. The state is measured before and after a fixed 
waiting period At to determine whether a decay of the 
excited state has happened. 



This three-step cycle is repeated typically several thou- 
sand times. The decay probability p is then determined 
as the ratio of D-state results in the second and the 
first state detections. For the calculation of the D5/2- 
statc lifetime T(^5/2) we use an exponential fit function 
(1 — p) = cxp{— Ai/T(5/2)}- For At we use the time in- 
terval between the ends of the two detection periods. 

Fig.Elshows the measured D5/2-level excitation proba- 
bility (1 —p) after several delay times ranging from 25 ms 
up to 5 s. A weighted least squares fit to the data yields 
the lifetime T(^5/2) = 1168(9) ms using the fitting func- 
tion described above, where the only fitting parameter is 
''"(5 /2) • The statistical error (in brackets) is the la stan- 
dard deviation. The fit yields xl — 0.47, indicating that 
the experimental decay is consistent with the expected 
exponential decay behavior. The least-squares method 
is justified by the normal distribution of the mean decay 
probabilities which is a result of the 'central limit theo- 
rem' of statistics. This was also verified using simulated 
data sets (see next section). 
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FIG. 5: D5/2-level excitation probability after delay times 
from 25 ms to 5 s plotted on a logarithmic scale. The solid line 
is a least squares fit to the data yielding ^5/2) ~ 1168(9) ms. 
The residuals (difference of data points and fit curve) of the 
fit are shown in the lower diagram. 



B. Systematic errors 

There are several types of systematic errors that may 
affect the lifetime result. In UHV single ion experiments 
the biggest error source is usually radiation which irradi- 
ates the ion due to insufficient shielding of room light 
or insufficient shut-off of laser beams. In our experi- 
ment, the strongest influence stems from residual light 
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at 854 nm. The influence of this radiation on the D5/2 
-level lifetime has been investigated extensively in . If 
radiation at 854 nm is present during the delay interval it 
may de-excite the D5/2-level to the ground state via the 
strongly coupled P3/2-level. This additional "decay chan- 
nel" artificially shortens the observed lifetime. The ob- 
vious source for residual 854 nm radiation is the 854 nm 
diode laser itself. In our experiment, it is eliminated by 
a fast mechanical shutter [33 which is closed during the 
delay interval. The 40 dB attenuation of the double-pass 
AOM which usually switches the 854 nm light was shown 
to be insufficient: In an earlier experiment without the 
shutter the lifetime was determined to 1011(6) ms jSj. In 
addition, the results observed without shutter were found 
to fluctuate by approximately ±50 ms depending on the 
specific AOM and diode laser adjustments. 

Another source for 854 nm radiation is background flu- 
orescence at 854 nm from the 866 nm diode laser. To 
eliminate this radiation an AOM in single pass with an 
attenuation of more than 20 dB was used to shut the 
866 nm beam. As the systematic lifetime error with- 
out AOM was found to be of the order of a few percent, 
this attenuation is sufficient. Note that this source of 
error cannot, in principle, be directly eliminated in the 
quantum-jump technique where 866 nm light must be ra- 
diated onto the ion continuously. In that case, the only 
way to correct for this systematic error is to measure at 
different light powers and extrapolate the lifetime to zero 
power which in turn implies a larger statistical error. In 
summary, radiation at 854 nm did not influence the mea- 
sured D5 /2 -level lifetime at the given level of statistical 
uncertainty. 

The D5/2 -level lifetime could in principle be also re- 
duced by transitions between the D-levels, i.e. by a Ml- 
transition stimulated by thermal radiation. The corre- 
sponding transition rate is given by W12 — Bi2p{v) with 
the Einstein coefficient for stimulated emission B12 and 
the energy density per unit frequency interval for thermal 
radiation p{v). With the rate of spontaneous emission 
^12 = {Sirhv^ / c'^)Bi2: W12 is rewritten as: 



With THz and A12 = 2.45 x 10"^ s"^ taken from 

[ssf we get 14^12 = 7.23 x 10^^ s^^ at room temperature 
which changes the D5/2 -level lifetime by much less than 
the statistical error of our measurement. 

Residual radiation could also induce lifetime- 
enhancing systematic effects. Both radiation at 393 nm 
(roomlight) or 729 nm (Ti:Sa laser, double-pass AOM 
attenuation of « 40 dB) can lead to re-shelving of the 
ion. This effect, however, leads also to a different decay 
function. It is modeled by a simple rate equation for the 
excited state population 

PD ^ -rpD + R{i - pd) (2) 

where F = l/r denotes the natural decay rate and R the 
reshelving rate induced by laser radiation. The solution 



of eq. (0) is of the form 

PD = Ae-^'* + B (3) 

with A = 1 — R/r', an offset B = R/V and the new decay 
rate T' — V + R. Thus, an offset B is the signature of a 
re-shelving rate. The result from fitting the experimental 
data in Fig. with the modified exponential fit function 
P is Tfit = 1165(10) ms and i?flt = 3(2) x 10"^ s'^. 

To evaluate the systematic error due to a repumping 
rate R we use the following technique: We generate simu- 
lated data sets from random numbers by considering the 
fact that the decay probability for a given waiting time 
is distributed binomially around a mean that is given by 
an exponential function with an expected mean decay 
time T(5 /2) . For these data sets we also take into account 
the particular experimental waiting times and number of 
measurements. In this way 'ideal' data sets are created 
that contain a purely statistical variation of data accord- 
ing to the experimental settings and that are free of any 
systematic errors. First, a fit of Eq. |(2Jl to such an ideal 
simulated data set yields an additional repumping rate 
R — with a standard deviation of AR = 3 x 10~^ s~^. 
Thus, the above fitted repumping rate i?fit is consistent 
with zero and not sufficiently significant to allow any con- 
clusion about the actual rate or the model, i.e. the sta- 
tistical error is too large for this small systematic error 
to be resolved in a fit to the data. Second, to obtain an 
upper limit for the systematic error of the lifetime due to 
a possible re-pumping rate we assume that such a rate 
-Rsim exists in the experiment. We then simulate data 
sets including the rate i?sim and the lifetime T(5/2) and 
fit these data sets with a normal exponential fit func- 
tion {1 — p) = exp (— At/Tsim). The deviation of the fit 
result Tsim from T(5/2) used for the simulation gives the 
systematic lifetime error. For i?sim = 3 x 10^"^ s^^ this 
systematic error is Ar = —3 ms. 

Another systematic effect that implies a different fit 
model is the state detection error. Even though the de- 
tection efficiency is close to unity, Poissonian noise in the 
PMT counts and the possibility of a spontaneous decay 
during the detection period produce a small error |39l |. 
The first error, i.e. the probability Si to measure the 
ion in the wrong state due to noise of the count rate, 
depends on the discrimination between S and D state in 
the electron shelving technique. Such discrimination is 
achieved by comparing the fluorescence count rate dur- 
ing the detection window with a threshold value. Proper 
choice of this threshold value leads to an error as small as 
ei = 10~^ which can be neglected for the following analy- 
sis. The second error, i.e. the probability S2 for a wrong 
state measurement due to spontaneous decay, also de- 
pends on the length of the detection window, fluorescent 
count rates for the ion in S and D states and the thresh- 
old setting. For the chosen parameters in the experiment 
we evaluate £2 = 10"'^. The measured excitation proba- 
bility is then (l-p),„cas (l-e2)(l-p)reai and implies a 
model function of the form (l—p) = (1 — £2)6"^^'. A flt 
to simulated data as described above yields a statistically 
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consistent limit for this detection error of £2 = 1 x 10~ . 
Again, it cannot be resolved by a fit to the measured 
data. From simulated data including an assumed detec- 
tion error of £2 = 1 x 10"'^ we get an upper limit of the 
systematic lifetime error Ar — 7 ms. 

In addition to radiative effects, non-radiative lifetime 
reduction mechanisms exist, namely, inelastic collisions 
with neutral atoms or molecules from the background 
gas. Two relevant types of collisions can be distinguished: 
Quenching and j-mixing collisions. Quenching collisions 
cause direct deshelving of the ion into the ground state. 
In the presence of high quenching rates lifetime measure- 
ments have to be done at different pressures. An extrap- 
olation to zero pressure then yields the natural lifetime. 
Measurements of collisional deshelving rates for differ- 
ent atomic and molecular species have been performed 
in early experiments 0, 0|. Ref. 0| finds specific 
quenching rates for Ca"*" of = 37 x 10""'^^ cm^s~^ 
for H2, and V\ = 170 x lO-^^ cm^s^^ for N2. Colli- 
sions may also induce a change of the atomic polariza- 
tion, a process called j-mixing or fine structure mixing 
where a transition from the D5/2 to the D3/2 state or 
vice versa is induced. These rates have also been de- 
termined in Ref. [3 to r]^ = 3 x 10-^° cm^s"^ for H2 
and rj^ = 13 X 10"^° cm^s^^ for N2. Such coUisional 
effects cannot be distinguished from a natural decay pro- 
cess. Collisional effects are most prominent in experi- 
ments with large clouds of ions or at higher background 
pressure. To give an upper limit of the effect in our ex- 
periment estimates of the constituents of the background 
gas must be made. If a backgrour id g as composition of 
50% N2 and 50% II2 is assumed [4^ and the pressure 
p < 2 X 10~^^ mbar in the linear ion trap set up is taken 
into account, an upper limit for the additional collision 
induced rate of 7 = 3 x 10~* is calculated. This effect 
is well below a 10^'^ relative error and can be neglected 
here. 

In summary, the result for the lifetime of the /2 level 
can be quoted as: T(^5/2) = 1168 ms ±9 ms (statistical) 
-3 ms (repumping rate) +7 ms (detection). 

IV. MEASUREMENT OF THE D3/2 STATE 
LIFETIME 

A. Measurement procedure and results 

For the measurement of the D3/2-level lifetime some al- 
terations in the pulse sequence are required (see Fig. |HJ). 
To populate the D3/2 state we use indirect shelving by 
driving the S-P transition at 397 nm and taking advan- 
tage of the 1:16 branching ratio into the D3/2 level. After 
a few microseconds the D3/2 level is populated with unity 
probability. 

Furthermore, because that level is part of the closed 3- 
level fluorescence cycle used for state detection its popu- 
lation cannot be probed with the state detection scheme 
described in the previous paragraph. In that sense the 



D3/2 level is not a shelved state. The method used here 
is that prior to state detection the decayed population 
is transferred to the D5/2 shelving state. The measured 
excitation probability of the D5/2 state divided by the 
shelving probability then corresponds to the decay prob- 
ability from the D3/2 level and the further analysis is 
analogous to the one in Sec. IIIII Shelving in the D5/2 
state is achieved by coherent excitation (7r-pulse). How- 
ever, it must be taken into account that the D3/2 state 
may decay into both Zeeman sublevels of the S1/2 ground 
state. Hence, two tt pulses from both sublevels are re- 
quired to transfer all decayed population to the D5/2- 
state. In our experiment we chose the two Aruj — 2 
transitions {rrij = —1/2 to rrij — —5/2 and rrij — 1/2 
to rrij = 5/2). The combined transfer efficiency P^r of 
the two pulses is determined in the first part of the pulse 
sequence (cf. Fig. El): after Doppler cooling the ion is not 
optically pumped into the Si/2('7ij = —1/2) ground state 
as usual but might populate both Zeeman sublevels. The 
measured transfer efficiency P.^ is used for calculation of 
the decay probability. 
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FIG. 6: Simplified pulse scheme for the D3/2 lifetime mea- 
surement. It consists of a measurement of the 7r-pulse trans- 
fer efficiency on the S1/2-D5/2 transition (prep, tt and detl); 
D3/2-state preparation (prep, s); waiting period At and state 
detection (tt and det2). For details of the pulse sequence see 
text. The waiting time is varied between 25 ms and 5 s. 

The measurement of the D3/2-level lifetime T(3/2) thus 
consists of a repetition of the following laser pulse se- 
quences applied to the ion. The sequence generally is 
composed of three steps (cf. Fig. E)) : 

1. Measurement of transfer efficiency i-^: state prepa- 
ration and Doppler cooling, consisting of 2 ms of Doppler 
cooling (397 nm and 866 nm light), repumping from 
the D5/2 level (854 nm light) and spontaneous decay 
into the Si/2(m=-l/2) or (m=-fl/2) Zeeman sublevel; 
TT-pulses on the S1/2 to D5/2 transitions (rrij — —1/2 to 
rrij = —5/2 and rrij = 1/2 to rrij = 5/2); state detection 
for 3.5 ms by recording the fluorescence on the S1/2 - 
Pi/2 transition with a photomultiplier. 

2. State preparation and shelving in the D3/2-level: 
2 ms of Doppler cooling (397 nm and 866 nm light), 
repumping from the D5/2 level (854 nm light) and op- 
tical pumping into the Si/2(m=-l/2) Zeeman sublevel 
(397 nm polarized light); shelving pulse at 397 nm 
for a few fj,s. 

3. Measurement of decay probability: free decay for a 
variable delay time; 7r-pulses on the S1/2 to D5/2 tran- 
sitions {rrij = —1/2 to rrij = —5/2 and rrij — 1/2 to 
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rrij = 5/2); state detection for 3.5 ms by recording the 
fluorescence on the S1/2 - P1/2 transition with a photo- 
multipher. 

The measured D3/2-level excitation probabiUty (I —p) 
is plotted vs. various delay times between 25 ms and 
2 s in Fig. [3 Again, the data have been fitted using 
the least squares method and the fit function 1 — p = 
exp(— Ai/r(3/2)). Here, p denotes the corrected decay 
probability p = Pex/P-m i-e. the detected excitation of 
the D5/2 level P^x corrected for the near- unity shelving 
probability (which is typically 0.98-0.99 on average). 
Since there is no correlation between and P^x in one 
experiment it is more appropriate to use for the correc- 
tion the average of Pn- for each delay time. The output 
from the fit is T(^z/2) = 1176(11) ms with — 0.68 indi- 
cating good agreement of data and exponential model. 
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FIG. 7: 733/2-level excitation probability for delay times from 
25 ms to 2 s plotted on a logarithmic scale. The solid line is 
a least squares fit to the data yielding T(3/2) = 1176(11) ms. 
The residuals (difference of data points and fit curve) of the 
fit are shown in the lower diagram. 



B. Systematic errors 

Also for this experiment systematic errors due to resid- 
ual light have to be investigated. The measured lifetime 
might be reduced by residual light at 866 nm or 850 nm 
present during the delay interval. This light would de- 
excite the D3/2 level via the P1/2- or P3/2-levels, respec- 
tively, and results in a faster effective decay rate. The 
main source of light at 866 nm is the corresponding diode 
laser itself which is switched with a single pass AOM with 
an attenuation of 20 dB. As this attenuation was found 
to be insufficient a fast mechanical shutter (cf. Sec. 
was installed which remained closed throughout the en- 



tire waiting period. The fluorescence background of the 
854 nm diode laser at 866 nm was found to be negligible. 
Light at 850 nm could de-shelve the ion via the P3/2 state 
and is expected to mainly originate from the fluorescence 
background of the 854 nm diode laser. For this laser, a 
double-pass AOM attenuation of 40 dB was proven to 
be sufficient since no effect on the lifetime could be mea- 
sured within a 5% error even if the laser was switched on 
at full intensity during the whole waiting time. 

Lifetime reducing effects are not obviously detectable 
because they only increase the decay rate while the func- 
tional shape of the decay curve remains the same, i.e. no 
offset is introduced for long delay times. The main con- 
cern in our experiment is the 866 nm light and extreme 
care has been taken to ensure that the shutter was indeed 
closed during the delay time. Before the 397 nm shelving 
pulse and between the tt pulse and the second detection a 
1 ms period has been inserted to allow for shutting time 
and jitter. During the lifetime measurements the correct 
shutting was checked by monitoring the transmission on 
photodiodes. In fact the shutters close fast within about 
400 /is but the start time is not well defined and jitters 
by about 500 /is. 

Lifetime prolonging effects can be induced by residual 
light at 397 nm or 729 nm which might re-excite the ion 
after it has already decayed. This re-shelving process 
can be detected as an offset as already pointed out in 
Sec, mil The 397 nm laser light is switched by two single 
pass AOMs in series (one before a fiber, one behind it, 
combined attenuation « 55 dB). Nonetheless, we used a 
shutter to exclude the influence of 397 nm laser light to 
the largest possible extent. To give a limit on the sys- 
tematic effect of any re-pumping source the same method 
as in section IIII Bl is applied. The experimental data is 
fitted with the rate model function Eq. |(2Jl yielding a rate 
of Pfit = 3(10) X 10^3 g~i^ xhe standard deviation of R 
for an simulated ideal data set is AP = 1.5 x 10"^ s~^ 
(with mean P = 0) , so again the rate is concealed by the 
statistical error. From simulations an upper limit for the 
systematic lifetime error of Ar = —2 ms is obtained. 

Another source of systematic error is vibrational heat- 
ing of the ion during the delay time. If, due to heating, 
the transfer efficiency P,r(Ai) is smaller after the waiting 
time than P,r(0) determined in the first part of the pulse 
sequence the correction for the transfer efficiency is too 
small and the actual decay rate is higher than measured. 
A TT pulse only has high transfer efficiency if the ion is 
in the Lamb-Dicke regime jMl^i i-e- Tj^n « 1 where 
is the Lamb-Dicke parameter and n is the mean phonon 
number. If the factor rj'^n becomes significant both the 
Rabi frequency Qn and the maximum transfer efficiency 
decrease as iln = i^o(l — v'^^), where fla is the coupling 
strength on the S-D transition. Taking the mean phonon 
number after Doppler cooling of n « 10 and the measured 
heating rate in the linear ion trap of dn/dt « 10 s~^ [s^ 
we can estimate the transfer efficiency after a waiting 
time of 2 s, Pt(2 s) = 0.98 if the 7r-pulse time was initially 
chosen to fulfill P7r(0) = 1. We experimentally checked 
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the degradation of transfer efficiency with waiting time 
by introducing a delay time At between the Doppler cool- 
ing pulse and the vr-pulse in the first step of the pulse 
sequence and subsequently performing a state detection 
measurement. Figure |S1 shows an average of various mea- 
surements of TT-pulse transfer efficiency P7r(Ai) vs. delay 
time At. A linear fit P7r(At) = 1 — a At to the data 



represented as an expansion 




FIG. 8: Average transfer efficiency of the 7r-pulses on the S1/2 
to D5/2 transition after various delay times between Doppler 
coohng and vr-pulses. 



yields a = -4(2) x 10"^ g- 



With simulated data sets 



including such a decreasing transfer efficiency Pn-(At) we 
determine a systematic lifetime error of At = —7 ms. 

Finally, the detection error is considered, in analogy to 
section IlII Bl From simulated data sets with a detection 
error of £2 = f x 10"'^ a systematic error for the life- 
time of At = +8 ms is found. Systematic errors due to 
collisional effects (quenching and j-mixing) can be again 
neglected as argued above. 

Summarizing the analysis, the lifetime for the 
level is given as: T(3/2) — 1176 ms ±11 ms (statistical) 
-2 ms (repumping rate) -7 ms (heating) -1-8 ms (detection 
error) . 



AB INITIO ALL ORDER CALCULATION OF 
THE D-STATES LIFETIMES 



We conducted the calculation of the 3d 



3/2 



4s 



^Si/2 and 3d ^05/2 - 4s ^Si/2 electric-quadrupole ma- 
trix elements in Ca+ using a relativistic all-order method 
which sums infinite sets of many-body perturbation the- 
ory terms. These matrix elements are used to evaluate 
the 3c?-level radiative lifetimes and their ratio. 

In this particular implementation of the all-order 
method, the wave function of the valence electron v is 



m^v mna 



iniirab 



I*.), (4) 



where $^ is the lowest-order atomic wave function, which 
is taken to be the frozen- core Dirac-Hartree-Fock (DHF) 
wave function of a state v. This lowest-order atomic wave 
function can be written as |<I>„) = aj,|0c), where |0c) rep- 
resents DHF wave function of a closed core. In equation 
a| and are creation and annihilation operators, 
respectively. The indices m, n, and r designate excited 
states and indices a and b designate core states. The first 
two lines of Eq. Q represent the single and double exci- 
tation terms. The restriction of the wave function to the 
first five terms of Eq. Q represents a single-double (SD) 
approximation. The last term of Eq. Q represents a class 
of the triple excitations and is included in the calculation 
partially as described in Ref. . We carried out the all- 
order calculation with and without the partial addition 
of the triple term; the results of those two calculations 
are labeled SD (single-double) and SDpT (single-double 
partial triple) data in the text and tables below. The ex- 
citation coefficients pma, Pmv, Pmnab, and pmnva are ob- 
tained as the iterative solutions of the all-order equations 
in the finite basis set. The basis set, used in the present 
calculation, consists of the single-particle states, which 
are linear combinations of B-splines 42]. The single- 
particle orbitals are defined on a non-linear grid and are 
constrained to a spherical cavity. The cavity radius is 
chosen to accommodate the 4s and 3c? orbitals. 

The matrix element of the operator Z for the tran- 
sition between the states w and v is obtained from the 
expansion using 



(5) 



The resulting expression for the numerator of the Eq. (0) 
consists of terms that are linear or quadratic functions 
of the excitation coefficients. We refer the reader to 
Refs. 10, 0, ^3 for further description of the all-order 
method. 

The numerical implementation of the all-order method 
requires to carry out the sums over the entire basis set. 
We truncate those sums at some value of the orbital an- 



gular momentum / 



max ^ ''iriax 



6 in the current all-order 
calculation. The contributions of the excited states with 
higher values of I which are small but significant for the 
considered transitions, are evaluated in the third-order 
many-body perturbation theory (MBPT). To evaluate 
those contributions, we carried out a third-order MBPT 
calculation with the same basis set and Imax, used the 
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TABLE I: Electric-quadrupole reduced matrix elements E2 
in Ca+ (in a.u.) calculated using different approximations: 
Dirac-Hartree-Fock (DHF), third-order many-body pertur- 
bation theory (Third order), single-double all-order method 
(SD), and single-double all-order method including partial 
triple excitation contributions (SDpT). The all-order data cal- 
culated with /max = 6 are listed separately. The contribution 
of basis set states with orbital angular momentum / = 7 — 10 
calculated using third-order MBPT is listed in rows labeled 
"Extr.". This correction is added to obtain the values listed 
in rows labeled "SD" and "SDpT" . 







Method 




oQ 


2 Pi /La 2q 
U3/2 - 4S Si/2 










Third order 


7.364 








7.788 






SDpTi„„^=6 


7.971 






Extr." 


-0.038 






SD 


7.750 






SDpT 


7.934 


3d 


2D5/2 - 4s 2Si/2 


DHF 


11.978 






Third order 


9.046 






SD,„„,=6 


9.562 






SDpT,„„^=6 


9.786 






Extr." 


-0.046 






SD 


9.516 






SDpT 


9.740 



"This value is the difference of the third-order result obtained 
with the same basis set as the all-order calculation (the number 
of splines N = 35/40 and Imax = 6) and third-order result with 
N = 70 and l„,ax = 10. 



all-order calculation and then repeated the third-order 
calculation with larger basis set containing the orbitals 
with I up to Imax = 10. The difference between these 
two results is added to the ab initio all-order results. 
The convergence of the MBPT terms with I is rather 
rapid; the differences between the third-order calcula- 
tions with Imax = 4,6,8,10 are 1.8%, 0.4%, and 0.1%, 
respectively. The last number is well below the expected 
uncertainty of the current calculation. Thus, the contri- 
bution from orbitals with Imax > 10 can be omitted at 
the present level of accuracy. The contribution from the 
excited states with Imax > 6 relative to the total value 
of the matrix elements is significantly larger for 4s — 3d 
electric-quadrupole transitions (about 0.5%) than for the 
primary ns — np electric-dipole transitions in alkali-metal 
atoms. We note that while the all-order matrix elements 
contain the entire third-order perturbation theory con- 
tribution there is no straightforward and simple way to 
directly separate it out (see Ref. ji^ for the all-order vs. 
perturbation theory term correspondence issue). Thus, 
we have conducted a separate third-order calculation fol- 
lowing Ref. The results of the third-order and the 



all-order calculations (with and without partial inclusion 
of the triple excitations) are listed in TableQ] The contri- 
bution from the excited states with orbital angular mo- 
mentum I > 6 calculated as described above is listed in 
the row labeled "Extr.". The all-order values corrected 
for the truncation of the higher partial waves are listed 
in rows labeled "SD" and "SDpT" . 

We also investigated the effect of the Breit interaction 
to the values of the electric-quadrupole matrix elements. 
The Breit interaction arises from the exchange of a vir- 
tual photon between atomic electrons. The static Breit 
interaction can be described by the operator 

Bij = — - a,-cxj + ^[a,-aj- {cXi ■ f^) {a^ ■ )] 

(6) 

where the first part results from instantaneous magnetic 
interaction between Dirac currents and the second part 
is the retardation correction to the electric interaction 
|47l |. In Eq. ©, a; are Dirac matrices. The complete 
expression for the Breit matrix elements is given in 4^. 
To calculate the correction to the matrix elements arising 
from the Breit interaction, we modified the generation of 
the B-spline basis set to intrinsically include the Breit 
interaction on the same footing as the Coulomb inter- 
action and repeated the third-order calculation with the 
modified basis set. The difference between the new val- 
ues and the original third-order calculation (conducted 
with otherwise identical basis set parameters) is taken 
to be the correction due to Breit interaction. We give 
the breakdown of the third-order calculation with and 
without inclusion of the Breit interaction in Table ^] 
The Dirac-Hartree-Fock values are given in column DHF. 
The random-phase approximation (RPA) values, iterated 
to all orders, are listed in column RPA. The third-order 
Brueckner-orbital, structure radiation and normalization 
terms are listed in the columns BO, SR, and Norm, re- 
spectively. The breakdown of the third-order calcula- 
tion to RPA, BO, structure radiation and normalization 
terms follows that of Ref. j4gj . The reader is referred to 
Ref. [4^ for the detailed description of the third-order 
MBPT method and the formulas for all of the terms. We 
find the Breit correction to the DHF contribution to be 
dominant, with the contributions from all other terms be- 
ing insignificant. The total Breit correction is very small 
and below the estimated uncertainty of our theoretical 
values discussed below. However, the Breit contribution 
to the ratio of the matrix elements is found to be small 
but significant owing to higher accuracy of the ratio. 

The procedure described above does not include a 
class of the Breit correction contributions referred to in 
Ref. as two-body Breit contribution 50]. To con- 
duct the study of the possible size of the two-body Breit 
contribution we calculated the Breit contribution to 10 
different electric-dipole matrix elements (6s — 6p, 6s — 7p, 
7s — 7p, 7s — 6p, and 5^3/2 — Qp) in Cs using the method 
described above and compared those values with the re- 
sults from Cs is chosen as "model" atom as it is a 
similar system compared to Ca+. The Breit contribution 
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TABLE II: The Breit correction to the third-order values of the 4s - 3d ^03/2 and 4s ^Si/2 - 3d ^05/2 electric-quadrupole 

matrix elements. The Dirac-Hartree-Fock values are given in column DHF. The random-phase approximation (RPA) values, 
iterated to all orders, are listed in column RPA. The third-order Brueckner-orbital, structure radiation, and normalization 
terms are listed in the columns BO, SR, and Norm, respectively. 





Transition 




DHF 


RPA 


BO 


SR 


Norm 


Total 


4s 


^Si/2 - 3d ^D3/2 


no Breit 


9.7673 


-0.0553 


-2.2136 


0.0621 


-0.1588 


7.4018 






with Breit 


9.7611 


-0.0552 


-2.2131 


0.0621 


-0.1589 


7.3961 






Difference 


-0.0062 


0.0001 


0.0005 


0.0000 


-0.0001 


-0.0057 


4s 


^Si/2 - 3d ^D5/2 


no Breit 


11.9782 


-0.0662 


-2.7006 


0.0756 


-0.1945 


9.0926 






with Breit 


11.9672 


-0.0662 


-2.7001 


0.0756 


-0.1946 


9.0820 






Difference 


-0.0110 


0.0000 


0.0005 


0.0000 


-0.0001 


-0.0106 



to Cs properties was studied in detail owing to its impor- 
tance for the interpretation of Cs parity nonconservation 
experiments. In Ref. |49|| . both one-body and dominant 
two-body Breit contributions have been taken into ac- 
count. We find the largest difference between our data 
and that of (4^ to be 25%. For most of the transitions, 
we either agree to all the digits quoted in 49] or differ by 
10% or less. Therefore, the two-body contribution was 
not significant for any of the Cs electric-dipole transitions 
that we could compare with. We agree with the values of 
the Breit correction to the DHF matrix elements given in 
Ref. exactly, as expected, since the two-body Breit 
contribution only affects the correlation part of the cal- 
culation. Thus, we found no evidence that the two-body 
Breit correction may exceed the already calculated one- 
body correction, especially considering the fact that the 
Breit correction to the lowest-order DHF value dominates 
the one-body Breit correction to the 4s — 3d matrix ele- 
ments in Ca^. Therefore, we assume that the two-body 
Breit contribution does not exceed the already calculated 
part. In summary, the omission of the two-body Breit 
interaction introduces an additional uncertainty in our 
calculation and we take the uncertainty to be equal to 
the value of the correction itself. Most likely, it is an 
overly pessimistic assumption based on the comparison 
with the calculation of the Breit correction to Cs electric- 
dipole matrix elements carried out in Ref. . 

Next, we use a semi-empirical scaling procedure to 
evaluate some classes of the correlation correction omit- 
ted by the current all-order calculation. The scaling pro- 
cedure is described in Refs. j44i l45l| . Briefly, the single- 
valence excitation coefficients pmv are multiplied by the 
ratio of the corresponding experimental and theoretical 
correlation energies and the calculation of the matrix el- 
ements is repeated with those modified excitation coeffi- 
cients. This procedure is especially suitable in this partic- 
ular study since the matrix element contribution contain- 
ing the excitation coefficients pmv overwhelmingly dom- 
inates the correlation correction for the considered here 
transitions. We conduct this scaling procedure for both 
SD and SDpT calculations; the scaling factors are differ- 
ent in these two cases as SD method underestimates and 
SDpT method overestimates the correlation energy. 



Table lnil contains the summary of the resulting matrix 
elements; the Breit correction is included in all values. 
We note that the scaled values only include DHF part of 
the Breit correction to avoid possible double counting of 
the terms (because of the use of the experimental corre- 
lation energy in the scaling procedure). The final values 
are taken to be scaled SD values based on the compar- 
isons of similar calculations in alkali-metal atoms with 
experiment 

EililEll- The uncertainty is calculated as 
the spread of the scaled values and ab initio SDpT val- 
ues. The uncertainty in the Breit interaction calculation 
is also included; it is negligible in comparison with the 
spread of the values. 

TABLE III: Electric-quadrupole reduced matrix elements E2 
in Ca^ (a.u.) 





Transition 


Method 


ab initio 


scaled 


3d 


2D3/2 - 4s 2Si/2 


SD 


7.744 


7.939 






SDpT 


7.928 


7.902 






Final 




7.939(37) 


3d 


^D5/2 - 4s ^Si/2 


SD 


9.505 


9.740 






SDpT 


9.729 


9.694 






Final 




9.740(47) 



We use our final theoretical results to calculate the life- 
times of the D3/2 and D5/2 states in Ca"*". The transition 
probabilities are calculated using the formula 



1.11995 X 10^8 |(^;||g||^,)|2 
A5 2j~n~ 



(7) 



where (uHQHw) is the reduced electric-quadrupole matrix 
element for the transition between states v and w and A is 
corresponding wavelength in A. The lifetime of the state 
V is calculated as 

r. = ^^^. (8) 

In both D3/2 and D5/2 lifetime calculations we consider 
a single transition contributing to each of the lifetimes. 
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The transition probabilities of other transitions (Ml D3/2 
- Si/2, Ml D5/2 - and E2 D5/2 - D3/2) have been 

estimated in Ref. [2a| and have been found to be 6 to 13 
orders of magnitude smaller that the transition probabil- 
ities of the D3/2 " Si/2 and D5/2 - S1/2 E2 transitions. 
Thus, we neglect these transitions in the present study. 
The experimental energy levels from Ref. [s^l are used 
in our calculation of the lifetimes. From the calculations 
we yield T(3/2) = 1196(11) ms for the D3/2-state and 
■''(5/2) — 1165(11) ms for the D5/2-state. These lifetime 
values are compared with experimental and other theo- 
retical results in Figs. El and 13 

The all-order calculation is in agreement with the 
present experimental values and recent experiments p3 . 
[23,123 within the uncertainty bounds. The present calcu- 
lation includes the correlation correction, which is large 
(23%) for the considered transitions, in the most com- 
plete way with comparison to all other previous calcula- 
tions m mm m 13 and is expected to give the most 
accurate result. It is also the only calculation which gives 
an estimate of the uncertainty of the theoretical values. 

In Ref. (2^, the issue of the theoretical ratio of the 
T(3/2)/r(5/2) lifetimes was raised. It appeared that there 
was a disagreement between previously calculated theo- 
retical ratios: Barton et al. I22|] quotes the values 1.0283 
1.0175 [23, and 1.0335l2l|. Such a disagreement 
appears to be rather puzzling since this particular ra- 
tio is far less sensitive to the correlation correction than 
the values of the corresponding matrix elements. Thus, 
we studied the value of the ratio and its uncertainty in 
detail. We list the values of the ratio of the D3/2 and 
D5/2 lifetimes calculated in various approximations in 
Table HVl The experimental energy levels from Ref. |5^ 
are used in all our calculations of the lifetimes for consis- 
tency. We include results with and without the addition 
of the Breit correction. As mentioned before, we find that 
the ratio does not change substantially with the addi- 
tion of the correlation correction; in fact, the correlation 
only contributes about 0.15% to the final value. Thus, 
we calculate the uncertainty in the ratio by considering 
the spread of the high-precision values of the ratio itself 
(SDsc, SDpT, and SDpTsc), rather than calculating the 
uncertainty in the ratio from the uncertainties in the in- 
dividual matrix elements. We also find that the while 
the Breit correction to the values of the matrix elements 
was insignificant at the current level of accuracy this is 
not the case for the ratio. In fact, the shift of the ratio 
values with the addition of the the Breit interaction is of 
the same order of magnitude as the spread of the high 
precision values as demonstrated in Table Hvl We take 
the SDsc value to be our final result for consistency with 
the calculation of the matrix elements. The uncertainty 
of the final value includes both the uncertainty in the 
correlation correction contribution and the uncertainty 
in the Breit interaction. As in the case of the individual 
matrix elements, the uncertainty in the Breit interaction 
is taken to be equal to the contribution itself. The Breit 
correction to the ratio is determined as the shift in the 



final ratio value due to addition of the Breit interaction. 



We compare our final theoretical value of the life- 
time ratio with the experiment and other theory in 
Tab le WJ The ratios of the other theoretical values 
El m nil are calculated from the numbers in 

the original publications; care was taken to keep the num- 
ber of digits in the ratio consistent with the number of 
digits in the values of the lifetimes or transitions rates 
quoted in the papers. First, we discuss the above men- 
tioned discrepancy of the theoretical ratios. Ref. lists 
the following ratios: 1.0283 "28^, 1.0175 '29^, and 1.0335 
[31.] . We have listed the data from the original publica- 
tions in Table which shows that the actual numbers 
with taking into account the number of digits quoted in 
the original papers should have been 1271/1236=1.028 
m, 1.16/1.14=1.02 m, and 1080/1045=1.033 |3l|. The 
first result is essentially a third-order relativistic many- 
body perturbation theory calculation with addition of 
semi-empirical scaling and omission of the some classes 
of small but significant third-order terms. It is very close 
to our third-order number 1.0286 from Table Hvl The 
next paper |29| quotes only 3 digits in the lifetime values 
(1.16s and 1.14s) so the accuracy is insufficient to obtain 
the fourth digit in the ratio. We note that the method de- 
scription in 29] is that of the non-rclativistic calculation 
and it is unclear how the separation to D3/2 and D5/2 life- 
times was made. The last calculation yields a larger ratio 
but that calculation has serious numerical issues such as 
taking only 20 out of 40 B-splines and including too few 
partial waves. It also omits all terms except Brueckner- 
orbital ones and possibly even third-order Brueckner or- 
bital contributions, which are large. The paper is not 
clear on the subject of the treatment of the higher-order 
contributions. Thus, we do not consider the result of ^| 
to be reliable. Therefore, there are essentially no incon- 
sistencies in the previously calculated theoretical ratios 
when the accuracy of the calculations is taken into ac- 
count. Our theoretical value of the lifetime ratio is higher 
than the experimental value. The spread of all values in 
Table including even lowest-order DHF values, is so 
small that it does not appear probable that any omitted 
Coulomb correlation or two-body Breit interaction can 
be responsible for the discrepancy. The only transition 
which can actually reduce the value of the theoretical 
ratio is the D3/2-S1/2 Ml transition. Thus, an accurate 
calculation of this transition rate will be useful in search 
for a theoretical explanation of the discrepancy. How- 
ever, the transition rate published in |38j| is extremely 
small {Ami — 7.39 x lO^^^s^^) and has to be incorrect 
by many orders of magnitude to affect the ratio at such a 
level which does not appear likely since the same calcu- 
lation gives a reasonably good (within 18%) number for 
the D3/2-S1/2 E2 transition rate. 
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TABLE IV: The ratio of the D3/2 and D5/2 lifetimes in Ca+ in various approximations. The lowest-order Dirac-Hartree-Fock 
results are labelled "DHF", third-order many-body perturbation theory results are in column labelled "Third", the results of 
the ab initio all-order calculation including single and double excitations are labeled "SD" , the results of the ab initio all-order 
calculation including single and double excitations with partial addition of the triple excitations are labeled "SDpT" , and the 
results of the corresponding scaled calculations are given in columns labeled "SDsc" and "SDpTsc", respectively. 





DHF 


Third 


SD 


SDpT 


SD,e 


SDpT,, 


Final 


No Breit 


1.0251 


1.0286 


1.0275 


1.0272 


1.0266 


1.0267 




With Breit 


1.0245 


1.0278 


1.0267 


1.0265 


1.0259 


1.0260 


1.0259(9) 



TABLE V: Comparison of the present values of the ratio of 
the D3/2 and D5/2 state lifetimes in Ca^ with other theory. 





Reference 


Value 


Theory 


m 


1.03 




[27j 


1.02 




[28J 


1.028 




[29j 


1.02 




[3^ 


1.033 




Present 


1.0259(9) 


Expt. 


Present 


1.0068(122) 



VI. DISCUSSION 

Figures |21 and 13 show an overview of the most recent 
experimental and theoretical results for the lifetime of the 
D5/2 and D3/2 states, respectively, in an chronological or- 
der. It is remarkable that the theoretical predictions scat- 
ter rather widely, with no noticeable convergence while 
the experimental results show a trend towards longer life- 
times in the recent years as more systematic errors are 
identified and stamped out. 

In comparison with previous work it can be concluded 
here that our lifetime result for the D5 /2 level agrees with 
and thereby confirms the most precise value of Barton 
et ai. We stress that this lifetime measurement is an 
independent check of earlier results as we used a different 
measurement technique. In addition, the result for the 
D3/2 level represents the first single ion measurement and 
reduces the statistical uncertainty of the previous values 
for the lifetime by a factor of four. 

For the calculated lifetimes we find excellent agree- 
ment of the theoretical all-order lifetimes with the ex- 
perimental results. Such agreement demonstrates the 
necessity of including partially the triple contributions 
to the all-order calculations for these types of transitions 
and confirms that scaling of the single-double all-order 
results, which is significantly simpler and less time con- 
suming calculation in comparison with ab initio inclusion 
of partial triple excitations, is adequate for these types 



of transitions. This is an important result for the evalu- 
ation of the accuracy of similar theoretical calculation in 
Ba^ which is important to parity violating experiments 
in heavy atoms. Such experiments are aimed at the tests 
of the Standard model of the electroweak interaction and 
at the study of the nuclear anapole moments. One of the 
features of most PNC studies in heavy atoms is the need 
for comparable accuracy of theoretical and experimental 
data. The current study is also of interest in regard to 
recently found discrepancy between the 5d lifetimes and 
the 6s — 6p Stark shifts in Cs Atomic properties of 
cesium were studied extensively by both experimental- 
ists and theorists owing to a high-precision measurement 
of parity non-conserving amplitude in this atom. Both 
of these quantities depend on the values of the 5d — 6p 
matrix elements. While those matrix elements are the 
electric-dipole ones rather than the electric-quadrupole 
ones studied here, the calculation itself as well as the 
breakdown of the correlation correction terms is very sim- 
ilar to the present calculation. Thus, the current study 
presents an important benchmark in the field of high- 
precision measurements and calculations. The study of 
the lifetime ratio demonstrated that the Breit interac- 
tion, which produces only a very small correction to the 
values of the actual matrix elements, is important in 
high-precision calculations of the corresponding matrix 
element ratios. 
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